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Abstract 

The Schmidt decomposition is an important tool in the study of quantum systems 
especially for the quantification of the entanglement of pure states. However, the 
Schmidt decomposition is only unique for bipartite pure states, and some multipar- 
tite pure states. Here a generalized Schmidt decomposition is given for states which 
are equivalent to depolarized pure states. Experimental methods for the identifi- 
cation of this class of mixed states are provided and some examples are discussed 
which show the utility of this description. A particularly interesting example pro- 
vides, for the first time, an interpretation of the number of negative eigenvalues of 
the density matrix. 
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1 Introduction 



Describing and quantifying entangled quantum states and preventing errors 
during quantum dynamical processes are important and, at this time, un- 
solved problems. Each of these has important implications for the develop- 
ment of reliable quantum information processing devices. In order to tackle 
these problems, it is common to broaden our knowledge and understanding by 
developing key examples. This is our approach here as we examine a particular 
class of quantum states. 
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This work was motivated by a desire to be able to identify and distinguish a 
certain class of mixed quantum states, and their properties, experimentally. 
This will rely, in part, on the existence of the Schmidt decomposition [1] which 
provides a canonical form for bipartite pure states. The Schmidt decomposi- 
tion is used to identify and quantify entanglement in bipartite quantum sys- 
tems [2]. Such systems are primitives for a host of quantum communication 
and computation protocols. However, such protocols are invariably subject 
to noise which diminishes their advantage over classical protocols. Noise, for 
most quantum systems, is difficult to identify and protect against, although 
there are several promising methods (see for e.g. [3] and references therein). 
Here we introduce a generalized Schmidt decomposition for a class of mixed 
quantum states which we hope will aide both with the problem of under- 
standing entanglement and our ability to correct for noisy quantum processes. 
Our decomposition does not retain all of the useful properties which make 
the pure-state version so important. However, it does allow us to devise some 
useful tools for measuring properties of an important class of states. 

The Schmidt decomposition is described by a set of real coefficients that is 
invariant under local unitary operations. All entanglement measures on pure 
states, such as the von Neumann entropy of a reduced density operator, can 
be computed from this set. However, this decomposition is known only to 
exist for general bipartite pure states (see for example [4]) and some multi- 
partite pure states [5,6]. Therefore, quantifying entanglement in terms of this 
decomposition does not work in general. For mixed states, several entangle- 
ment measures exist, most of which are difficult to calculate, though some 
interesting special cases for bipartite systems can be solved. For example, for 
two qubits one can calculate the Entanglement of Formation (EoF) [7] which 
is the amount of entanglement required to form a particular state. It is also 
known how to calculate the EoF for Werner states [8] , isotropic states [9] and 
rotationally invariant states [10]. However, at this time there is no canonical 
Schmidt decomposition for mixed states and no efficient method by which to 
analytically compute the entanglement of general mixed states. 

One might anticipate that a generalization of the Schmidt decomposition 
would aid in the description of entangled states. One such generalization is 
given by the Schmidt number [11], which is equal to the maximum Schmidt 
rank (or number of Schmidt coefficients) in a pure state decomposition of a 
mixed state, minimized over all decompositions. This quantity constitutes the 
minimum Schmidt rank of the pure states needed to construct a state, and is 
an entanglement monotone [11]. Here we consider another special case which 
is a Schmidt decomposition for depolarized pure states (DPS) which are those 
states obtained by mixing the identity operator on the state space with a sin- 
gle pure state. These have many interesting properties and have been studied 
in the literature since these states are fairly easy to manipulate. For example, 
one may compute properties such as channel capacities [12,13], entanglement 
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(specific instances) [14,15], and more recently, it has been shown that noisy 
operations may be turned into depolarizing operations [16]. The set of DPS 
which we define here includes, not only pure states which have undergone a 
depolarizing operation, but also states which, if initially decoupled from their 
environment, cannot be obtained in this way. All states in our DPS class can 
be brought into a similar canonical form using local unitary operations. 



The DPS are important to understand in part because they have a fairly simple 
form. This form has 2D — 1 real parameters as opposed to — 1 parameters 
for a generic mixed state in a D dimensional Hilbert space. They are also 
important to understand because any map can be brought to the depolarizing 
form by a simple sequence of quantum operations. Therefore a complicated 
quantum computing process in the presence of noise can be brought into this 
form which produces states with relatively few relevant parameters. This al- 
lows a direct comparison of inequivalent noise processes by projecting them 
into the same class. 



In this article we discuss methods for experimentally determining whether this 
form has indeed been produced. We find expressions for the fidelity and the 
trace distance for this class of mixed states, and are also able to show that 
the negativity is more easily quantified for bipartite DPS. More importantly 
perhaps, we provide a hound for the number of negative eigenvalues for bipar- 
tite DPS and show that the number of negative eigenvalues can indicate the 
type of entanglement present in the system, e.g. qubit-qubit vs. qutrit-qutrit. 
These results support a limited form of a conjecture by Han, et al. [17] about 
the maximum number of negative eigenvalues for a bipartite state. We em- 
phasize that our results provide an experimentally detectable qualitative and 
quantitative measure of entanglement. 



The paper is organized as follows. In Section 2.1 we review the coherence vector 
parameterization of the density operator. In Section 2.2 we provide a geomet- 
ric interpretation of DPS in terms of the coherence vector parameterization. 
Section 2.3 demonstrates that there exists a type of Schmidt decomposition 
for depolarized pure states when there exists a Schmidt decomposition for 
the corresponding pure state. In Section 3 we provide two ways in which to 
identify these states experimentally, and describe physical maps which give 
rise to DPS beginning in an unknown pure state. In Section 4 we discuss the 
insight that we gain into bipartite entanglement given our construction. We 
then conclude with a summary and some open questions in Section 5. Some 
examples of the formalism are given in Appendix A. 
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2 Schmidt form for DPS 



In this section we provide several forms for the DPS which will be used for 
various calculations in later sections. 

2. 1 The coherence, or Bloch, vector 

The generalized coherence vector, or Bloch vector representation [18,19,20,21] 
will provide a convenient geometric picture for several parts of our argument. 
For a two-state system the description is well-known. The general case for an 
D-dimensional system is presented here and the two-state system will be seen 
to be a special case. 

Any density operator p belonging to the set of bounded linear operators B{p) 
with Hilbert space dimension dim(7Y) = D, can be expanded in a basis con- 
sisting of the identity operator and an operator basis for su{D), the algebra of 
SU{D). Throughout this work, we represent the latter with a set of Hermitian, 
traceless matrices, {Xi}fLi^ which obey the following orthogonality condition 



The commutation and anticommutation relations for this set are summarized 
by the following product formula 



Here, is the D x D unit matrix, the Cijk are the structure constants of the 
Lie algebra represented by these matrices, and the dijk are referred to as the 
components of the totally symmetric "rf-tensor." 

The density matrix for an D-state system can now be written in the following 
form 



where n-\ = X]f ^ niXi. For D > 2 the following conditions characterize the 
set of all pure states. 



Tr(AiAj) = 26ij. 



(1) 



AjAj — —SijlD + icijkXk + dijkXk- 



(2) 




(3) 



n-n = l, and n-kn = n, 



(4) 
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where the "star" product is defined by 



[a -kb)k = ]l ^ _ ^ dijkaibj. (5) 



For D = 2, the condition n ■ n = 1 alone is sufficient [22]. Note that 



r^. = ^^(^ Tr(pAO. (6) 



To recover the case of the two-state Bloch sphere, note tliat the constants 1/D 



and y D{D — l)/2 reduce to 1/2 and 1 respectively, and the dtjk are identically 
zero, so the second condition in Eq.(4) is not required. In fact, as noted, it 
cannot be satisfied. 



2.2 Depolarized Pure States 



Throughout this paper we focus on a special class of mixed states, the depo- 
larized pure states (DPS). Such states are given by a (not necessarily convex) 
sum of the identity operator and a pure state: 

Pd = {l~p)^lD + PP, (7) 



for p' some pure state. By the unit trace and positivity conditions, we have 
—1/{D — 1) < p < 1. Letting cd = ^Jd{D — l)/2, we may rewrite this in a 
more suggestive form as 

Pd = j: (1d + CDpn-\) . (8) 



We note that for D > 2 the characterization is unique, i.e. pd corresponds 
to a depolarized form of a single pure state with coherence vector n. This 
is because the condition n -k n = n demands that both n and —n cannot 
correspond to physical pure states. Hence, any vector of the form pn has a 
unique purification, namely n. For D = 2 this is not the case because both 
n and —n correspond to pure states. From this latter form, we may interpret 
the DPS as arising from the affine map: n h-» pn, on the — 1 dimensional 
real vector space of coherence vectors. 

This provides a geometric description of the set of depolarized pure states. 
The space of DPS with a given p is isomorphic to the set of pure states (for 
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D > 2). (See for example [23] and references therein.) To see the geometry 
more exphcitly, note that the DPS can be written in the form 



Pd 



D 



pW 



V 



Note that the same matrix W will diagonalize both the pure state and the 
depolarized pure state. 

We will make use of this form to analytically compute the trace distance and 
fidelity between two DPS. The fidelity between two density matrices is defined 
by 



We consider two DPS both in a D dimensional Hilbert space, 



(9) 



where — ^ P < 1 and the overlap in their purifications is i^d^I/) , |$)) 
1$) p = /. The (square root) of the fidelity is 



F{pd,ad) = {D-2)y/E+Y.± 



2a+{b+2c)f+d+b{l~f) _ 
2 



where the parameters are given by: 



(10) 



a = 
b = 


(l-p)(l-g) 
D ' 








c = 




-l)p 


+ m-p)- 




d = 


[l-q+Dqf] 


{{D- 


2)p + 2-2y 


\{D-l)p+l){l 




+^(^ 


K{D- 


-l)p+l)(l- 


-p) - (1 -p)). 
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Fig. 1. Distance measures between two depolarized pure states (DPS) in a 
D = 9 dimensional Hilbert space with equal polarizations p = q. The Bu- 
res metric, trace distance, and fidelity satisfy the following inequalities [24]: 
B{p,a)'^/2 < D{p,a) < y/l — F{p, a) and surface plots of these three functions 
are shown. The value Pmin = —1/80 is the minimum polarization of depolarized 
pure states which are obtainable from a completely positive map. (We call these 
physically depolarized pure states, cf. Sec. 3.4.) 

The square root of the fidelity can be converted into a metric, specifically 
the Bures metric via B{p,a) = \ 2 — 2 JF(p, a), and an angle A{p,a) = 



cos~ y F{p, a). In the pure state case, the Bures metric is the Euclidean dis- 
tance between the two pure states with respect to the norm on the state 
space and the cosine of the angle between the states is the overlap. The Bures 
metric between two mixed states can be interpreted as the Euclidean distance 
between purifications of the mixed states minimized over all such purifications. 

One can also compute the distance (in the trace norm) between two mixed 
states. The distance is 



D(p, a) = -Tr|p — a 



tr; 



where the trace norm is defined lOLr = VOW. For the two DPS 



D(pd,ad) 



iD-2)\q-p\ 
D 



Ed 



iq~p]{l^D/2) 
D 



(12) 



The distance between two mixed states with the same coherence vector mag- 
nitude p = q is simply D{pd,ad) = (1 — f)\p\- The distance and fidelities of 
equally polarized pure states are plotted in Fig. 1. Notice that beginning in 
a pure state, i.e. p = 1, the distance and Bures metric between states with 
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/ < 1 will decrease under a depolarizing map until both states are mapped to 
the identity. For even stronger maps, i.e. p < the distance begins to increase 
again. As discussed in Sec. 3.4, the minimum value of polarization obtain- 
able by a physical map acting on input pure states is Pmin = — 1). 
At this value, the distance between the output states is D{pd,ad) = -^rzi- 
Thus we find that the distance (fidelity) between initially pure states is not a 
monotonically decreasing (increasing) function of the depolarization strength. 



2. 3 Schmidt Decomposition For A Pure Bipartite State 



To fix notation, let us recall the Schmidt decomposition for a pure state of a 
bipartite quantum system m D = xDb dimensions with subsystems A and 
B which have dimension and Db respectively. Without loss of generality, 
we will assume that Da < Db- Now, let 

p'AB = \^m, (13) 



where 

|vI/') = E«^mI'/'/)aIV'm%- (14) 

According to the Schmidt decomposition [1], there exist unitary matrices U 
which acts only on the first subsystem, and V which acts only on the second 
subsystem, such that can be written in the form: 

I*) = E^. (15) 



where the set {|0)^} ( fornis an orthonormal basis for A (B). In other 

words, there are local unitary transformations, U and V such that 

U®Vm = j:b,\<P,)A\^,)B, (16) 



where 

^i^jv = Uijaif^Vf,^. (17) 

and U, V can chosen so that the {bj} are real and positive. We will say that 
Qifj, is "diagonalized" [25] by the local unitary transformations U and V. The 
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reduced density matrices pA = TrB(p) = I't'j) a = TrA(p) 

b'j I'ipj)^ {"ipjl have the same eigenvalues b'j. 



Now, let us consider the density operator 



P 



D 



\ik,iil3 



1 ® 










« / 


we 


see that 


if the matrix 



:i8) 



Defining the matrix a^^ = q^, „v. ...^ ...^...^ 

diagonalized by U and V, then a^^ can be diagonalized by the same U 
and V. 

Pd = ^{^-p)+pT.hh' ® l^.-X^/l • (19) 

hi' 



Therefore, there exists a preferred local unitary basis for depolarized pure 
states and we refer to this preferred basis as the Schmidt decomposition for 
DPS. 

Furthermore, we can provide a relationship between the eigenvalues of the 
reduced density matrices for the two subsystems. Tracing over the subsystem 
B produces 

Pd^ = TTsiPd) = ^a-p)+pJ2b' • (20) 



Now, let us suppose that there are n non-zero eigenvalues of = Ti b{p') 
given by {b'j} with J27=i bf = ^- (Alternatively, we could let the sum go to Da 
noting that for some i, the eigenvalue could be zero.) Then the eigenvalues of 
are {^{^ —p) +pbi}i=o U {(1 — p)/Da}^^~'^- Tracing over the subsystem 
A produces 



TrA(p.) = ^(1 -P) +pE&? I^.X^.I • (21) 



The eigenvalues of p'^ are given by {-^^^(1 — p) +pb'f}^^Q L\ {{1 — p)/Db}^'^ 

There are two properties of the Schmidt decomposition which make it par- 
ticularly useful and are properties which one would want to preserve in any 
generalization. It specifies (i) preferred bases of (ii) bi-orthogonal states. It is 
clear that property (i) is retained for DPS. This relies on the fact that it is 
unique for pure states [4] barring a degeneracy in the spectrum of one of the 
subsystems. 
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The Schmidt decomposition for general bipartite DPS is the preferred basis 
which agrees with the pure state Schmidt decomposition counterpart of the 
DPS. This definition clearly retains the property (i) and it can be generalized 
to any system with a corresponding pure state Schmidt decomposition. For 
example those described by a multipartite Schmidt decomposition [5,26] will 
also have corresponding set of DPS. 

Can this preferred basis be used to quantify the entanglement of the system? 
Certainly this is not true for the entropy of the partial trace as can be seen by 
considering the extreme case where p = 0. However, we will discuss how the 
Schmidt form helps identify and distinguish certain types of entangled states 
in Section 4. 



3 Preparation and Identification of DPS 

It is now pertinent to ask, how does one know if a density matrix describes a 
system whose state is in the class DPS? Is there a way to characterize maps 
which give rise to these states? This section will provide the answers to these 
questions. 



3.1 State Tomography 

Using state tomography the elements of the density matrix may be determined. 
There are several ways in which to do this, some of which are more efficient 
than others. For our purposes, it is assumed that state tomography data has 
been collected and from it the coherence vector n determined, for example via 
Eq. (6). 

From Eq. (4) the coherence vector of a pure state satisfies n-k n = n. For a 
DPS, n — > pn, so that n-k n ^ p^n-k n, etc. From these relations, it is clear 
that all invariants described in [20] can be calculated by noting that for a 
DPS [fi*]''fi ■ n = p^^"^. Therefore the invariants reduce to the simplified form 
which is obtained by replacing n with p everywhere and neglecting the types 
of products. In other words, 

n ■ n = P^i 

n-kfi ■ n = P^i 
n-k n-k n ■ n = p^, and so on. 
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These conditions may be stated equivalently, and more succinctly, as 

n-'n = p'^, and n'kn = pn. (22) 



Note that, similar to the pure state conditions, these two conditions alone 
determine the set of eigenvalues for the density operator. 

Note also that the DPS with p < and with p > can be distinguished 
with the unitary invariant n-kn-n (provided D > 2). Hence given some prior 
certificate that the state is a DPS, we obtain complete spectral information 
from the measurement oi n ■ n and n-kn-n including the value of p. 

Alternatively, one may examine the eigenvalues of the system. If the eigen- 
values are given by a,b,b, ...,b and having a + {D — l)b = 1, then the system 
is in the class DPS. Notice that the spectrum of the bipartite density matrix 
can be used to define the class and this is unchanged by a global unitary 
transformation. 



3.2 Invariant Polynomials 

Another measurement process which will efficiently identify the DPS is due 
to Brun [27]. He showed that, in principle, the invariants Tr(p'") could be 
measured efficiently. From these, the eigenvalues may be determined. 

Let S be an operator which cyclicly permutes states of the system: 

s i^i) m ■ ■ ■ m = m 1^2) ■ ■ ■ i^n-i) , (23) 

then 

Tr(^p®'") = Tr(p'"). (24) 
To show this is quite straight-forward. Let 

i i 

be an orthogonal (Tr(p*/>') = 5'-') pure-state decomposition of the density 
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matrix. Then 



Sp»^ = Pi® P2®---Pm 



(26) 



Taking the trace simply produces a series of Kronecker deltas which force all 
Pi to have the same index so that 



(27) 



A physical implementation of this measurement can be realized using an in- 
terferometer type circuit. This works by preparing an ancilla qubit a in the 
state \+x)a, (|±x)a = 1/V^{\0)a ± |l)a)), and applying a sequence of m - 1 
controlled- SWAP gates between the ancilla and pairs of copies of p: 



U 



n io). 

j=0 



li...m + (1| ® SWAP(m - J, m - J - 1), 



where SWAP(r, s) = Y.^jZo \i)r 01 ® \ j)s (^1- Each controlled-SWAP gate can 
be implemented using 0{D'^) elementary two qudit gates [28]. A final measure- 
ment of the ancilla in the \^x)a basis gives measurement outcomes m = ±1 
with probability P{m = ±1) = Ml + Tr[p™]). 



Since the above result really only depends on the production of the appropriate 
delta functions, in practice, any cyclic permutation which is not the identity 
could be used. In fact, it need not be cyclic as long as there is no invariant 
subspace. 



One may suppose that a particular experiment may provide for a more efficient 
measurement using the polynomials. However, it may also be the case that 
some state tomography data is available or some partial information about 
the state is known. In either of these cases, it is relevant to note the Tr(p"*) 
and the coherence/Bloch vector are directly related [20,21]. 
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3.3 Efficient determination using local measurements 

Knowing that a system is in a DPS enables the determination of the eigen- 
values of pd with the determination of Tr(p^) and Tr(p^) alone. However, if 
we do not know whether or not the combined system is in a DPS, a natural 
question is, how could this be determined? Generically this could be achieved 
by measuring the full spectrum of the state as outlined above by performing 
D measurements over a total oi D{D + l)/2 identically prepared copies of the 
state. For bipartite systems, simpler measurements on the subsystems A and 
B can reveal partial information about the state. While such information is 
not sufficient to verify that the joint state is of DPS form, one can check for a 
violation of the consistency relations given in Sec. 2.3 that can rule out that 
possibility. For example, one can measure the spectrum of the reduced states 
PAiPb and verify that the two sets of eigenvalues are equal up to the scaling 
which depends on the dimension. Another, perhaps simpler, measurement is to 
verify that the density operators are full rank. If one reduced state was found 
to have rank less than its dimension, for example by obtaining a zero value in 
a projective measurement, then the corresponding combined state pab could 
not be a DPS. Furthermore, for Db > Da + 2, there must exist a degenerate 
subspace of the subsystem B of dimension Db — Da- If this is not present, the 
system cannot be in a DPS. 

3.4 Physical depolarization channels 

It is natural to ask if all states pd can be generated by beginning in a pure 
state p' and applying a physical map which depolarizes that state to the form 
pd- It turns out that this is not always possible. Rather, according to the value 
of p, there is a continuous subset of DPS that cannot be so generated. To see 
this, consider the class of maps 

Sp{p) = {l-p)^lD + PP = P- (28) 

In ref. [32] it was shown that maps Sp with — l/(-D — l)<p<l are positive, 
but only those with — 1/(D^ — 1) < P < 1 are completely positive. Completely 
positive maps (CPM) are those maps which act as the identity operator on 
an environment when the input is a tensor product state of the system and 
environment. Such maps are deemed to be physically allowed maps acting on a 
system which is uncorrelated with its environment. (However, some dynamics 
need not be completely positive [33,34,35].) The map £^p=_i/(i52_i) is termed 
the universal inverter as it outputs the positive operator closest to being an 
inversion of the coherence vector of an arbitrary input state. Given this de- 
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marcation we classify all states pd which are obtainable from a single copy 
of the (generically unknown) pure state p' via a CPM to be physically depo- 
larized pure states (PDPS). The criterion that the map act only on a single 
copy is emphasized because more powerful operations are possible using mul- 
tiple copies. For example, given an infinite number of copies of a pure state p' 
one CPM is to perform state tomography and from the classical information, 
synthesize pd exactly. 

One can synthesize any positive density operator p5 in a -D dimensional 
Hilbert space by preparing an entangled state of the system with a D di- 
mensional ancilla a and tracing over the ancilla. Namely, given an eigen- 
decomposition of the state ps = J2jPj \'^j)s (^il' prepares the pure state 
\'^)sa ~ \/Pj\4'j)s b)a' traces over the ancilla. Clearly this synthesizes 
any DPS. Yet, for an initially uncorrelated system and environment, the trans- 
formation is generically non-linear. Often it is the case that one is interested 
in generating a PDPS output given an unknown pure state p'g as input. This 
can be useful to drive noisy maps with many parameters on pure states, to a 
standard form of a quantum channel with only one parameter, namely p. We 
now discuss two protocols to do so. 

The first method is a variant of a construction in [32]. Here one performs joint 
operations on the system and two ancillary qudits ai and 02 each of dimension 
D. The initial state is a tensor product state of the system S and the ancillae: 

P = Ps ® \x)a,a, {X\ (29) 



where = a |0)„^ ^ Ef=V and |$+) = ^ Ef=V |j) |j) 

is the maximally entangled state. The parameter a can arbitrarily be chosen 
real. We are interested in the case where the system itself is composed of two 
parts A and B but for simplicity we treat it as a single system whose Hilbert 
space is spanned by the orthonormal states {\j)g}f=Q- The next step is to 
apply a unitary composed of pairwise coupling gates between qudits: 

Usa.a, = [nr^Xi® {jm?^'Xj,^' ® (/!] ^^^^ 

[Uf~'\j)s{j\^zim?~'\j)sij\®zQ. 



Here the unitary operators are defined X = Y,j \ j + 1) (j | and Z = Y,j e''^'^^^^ \ j) 
The action of this unitary on a pure state input for the system is Usaia2 li^) s l^)ai 
a I'i/')^! 1*^^)0102 1^ l^)ai \^^) Sa2- Upon tracing over the ancillae, the residual 
system state is then: 

Ps = T^T^aiailUsaiaiPUsaiail ^^-^^ 

= {i-\p?)p's + m'¥, 
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where by the normahzation constraint on the state < < D"^ jiTP' — 
1). Hence, by varying the parameter /3, one can reahze any PDPS. 

A second protocol for generating PDPS works by using stochastic unitaries 
to randomize a quantum operation E on an input state [16]. The degree to 
which the map S, acts trivially determines the depolarization parameter p and 
the randomization guarantees that the map takes all inputs p to the standard 
form Pp. Specifically, one randomly picks a unitary V G U{D) and applies U 
before and W after a trace preserving, CPM S on the state. The result is 

S'{p)=JdUUS{WpU)W 



where dU is the invariant Haar measure on U{D). Here < / < 1 quantifies 
the identity portion of the map, i.e. / = ($"'"| |$+) where is the Choi- 
Jamiolkowski representation [36,37] of the map 8. Such a representation arises 
by first writing a trace preserving CPM on B{l-is) in a particular operator- 
sum decomposition as £{p) = J2m7n,rn',n'=oEni,n-,m',n'X"-Z"'p{X"-' Z'^'V ■ The 

state E e B(Hs X Hs') given hj E = Em^|m',n'=0 -^^n-m'n' \^m,n)sS' {^m',n'\ 

expanded in the orthonormal basis {\^m,n)ss' ~ -^s^s" \^^)ss'}^ then the 
Jamiolkowski representation of S. This follows by virtue of the relation E^ = 

A simple way to generate a particular PDPS is as follows: 

• Begin with a pure state p'. 

• Pick a unitary U G U{d) at random and apply it to the state. 

• Apply a quantum operation with Jamiolkowski fidelity /; for example, the 
single qudit unitary V = e^"(^^+^A) which has / = ^^-j Efio"^ (,i2acos{2^j/DA) 
Another option is to apply the operator Xa with probability 1 — / and with 
probability / do nothing to the state. 

• Apply t/^ to the state. 

The resultant state is pd withp = jyTzr- practice, for the stochastic process, 
it is not necessary to pick a unitary uniformly at random, rather one can pick 
a random unitary from the finite set G = G \ 1, where G is the Clifford 
group. The latter is defined as the group which leaves the Pauli group P = 
|gi2^fc/Dj5^a^b. a,b,keZD} invariant under conjugation. 

We stress that both of the above protocols require performing entangling op- 
erations between the subsystems A and B. This is because in both cases, it 
is necessary to implement the Pauli operators Xs and Zs which cannot be 
written as local unitaries on A and B alone. This emphasizes the fact that the 
depolarizing map is a map on the joint space, it cannot be realized by sepa- 
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rately depolarizing each party. In fact the action of individual depolarization 
is a map with 4 real parameters: 



1 1,1 

"ooPab + "oiPa ® + "107:;- ® Pb + "117^ 

L>B L>A 



which is not the desired form. 



4 Entanglement of DPS 



Given the results of Section 3, we can determine experimentally whether the 
state has the form of a DPS or not. From this information we find the negative 
eigenvalues which provides a sufficient condition for the existence of entangle- 
ment in a mixed state. For a two qubit system, or a qubit-qutrit system the 
criterion is both sufficient and necessary. 



4-1 Partial Transpose 



Since partial transpose is independent of local unitary operations, we can 
compute it for the Schmidt form of a depolarized state. The explicit form of 
the partially transposed state is: 

= {l-p)lA^+pY.fS'b]\ct^,)M®\^3)B{i'3\ (33) 

Da-1 

i<i'=o 



where we introduced the orthonormal states: |±jj') = 
Notice that this form is diagonal. 



4 ■ 2 Negativity 



For states p G B{Ha ^ 'Hb) with Da{b) = dim7YA(B) the negativity N{p) is 
defined [29]: 
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where, again, without loss of generahty we assume Da < Db- The func- 
tion is real valued and normalized to lie in the range [0,1]. The argument 
is the partial transpose of p with respect to subsystem B, which in 
a coordinate representation with p = PijA'j' \i) a a {'i'l ® \j)BB{j'\, is 

p^'^ = Pij,i'j'\i) AA{i'\ ® IjObbO'I- While it's action is locally basis 

dependent, the eigenvalues of p"^^ are not, and the negativity counts a nor- 
malized sum of the norm of negative eigenvalues. Because any separable state 
can be written as a convex sum of products of partial density operators, and 
hence has eigenvalues invariant under partial transposition, negative eigenval- 
ues are a sufficient hut not necessary condition for the presence of bipartite 
entanglement in p. States with p^^ > but not separable are known as bound 
entangled states because that entanglement cannot be distilled. 

From Eq. 33 the negativity is quickly found to be: 



Af{pa 



Da-1 



Pbjbj'l . I B y-J-J 



^-pbjby\)-l 



l~p_ 
2Db' 



(35) 



All that is required for M{pd) > is that one of the terms inside the absolute 
value be negative or p > j-n^l^j^i for some pair of Schmidt coefficients bj,bj'. 

Notice, that since bjbjf < 1/2, then for p < jjj^:^, -f^ipd) = 0. It is also true 
that for p < jyj^^, the state is separable [30]. 

However, let us note that, from the diagonal form, we can extract more infor- 
mation. Any quantifier of entanglement, such as the EoF, or negativity, tells us 
only how entangled a state is. For quantum information purposes, we may like 
to know what type of entanglement is present in the system. For example, for 
distillation protocols, we may want to know if a type of qutrit entanglement 
is present. This is particularly relevant given that some quantum information 
protocols require entangled qudits. Let us consider what we may discern from 
Eq. (33). 



4-3 Number of Negative Eigenvalues 



The number of negative eigenvalues of the partially transposed joint state 
provides a sufficient condition for stratification of the pure state entanglement. 

Before addressing this point, recall from Sec. 3 that given some prior knowledge 
that a bipartite system is in a DPS, one may obtain the eigenvalues, i.e. the 
set {bi}, as well as p from the spectrum of one of the local density operators 
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alone, e.g. from pA = TTB{pd)- In what follows, it is assumed that the state is 
in a DPS and that p and {bi} have been determined. 



From Eq. (33), the eigenvalues of the partially transposed density operator 
will be 

{{(1 -p)i + pb',}fS\ {ii-p)j5+PWf<r\ ^3g) 
{{i-p)i-pb,b'^}f-f'} 



Note that the number of negative eigenvalues is bounded above by \^ . For 
two qubits this means that the maximum number of negative eigenvalues is 
one. For two qutrits, the maximum number of negative eigenvalues is three, 
etc. Note that for a maximally entangled state of two identical systems of 
dimension D^, 

1 Da-I 

= ^ E K^) ' (37) 

V^A i=o 



and symmetry requires that there are either Da negative eigenvalues or none. 
This result supports the conjecture by Han, et al. [17] that for the maxi- 
mum number of negative eigenvalues for a bipartite entangled mixed state is 
DAiDA - l)/2. (Recall Da<Db.) 

For example, consider = 9, and Da = 3 = Db- The eigenvalues of the 
partially transposed density operator are 

{^+Pbt)\^+pbt),(^+pbt), 

By inspection, any of the last three will be negative when 
1 

^ ^ %jbj, + 1' 



for a given j, j' as is consistent with the general requirement that the state 
be entangled according to the negativity. However, note that if cor- 
responds to a Bell state, then ^1 = ^ = ^2 and 63 = 0. This implies that 
there is at most one negative eigenvalue which occurs when p > 2/11. Now 
consider the maximally entangled two-qutrit state, 61 = 62 = ^3 = l/"\/3 (or 
any state locally equivalent to an SU(3) singlet). In this case, when p > 1/4, 
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all of the last three eigenvalues are negative. Clearly this cannot happen for 
a two qubit density operator since, at most, one eigenvalue is negative. 
The difference in the number of negative eigenvalues therefore provides a suffi- 
cient condition for distinguishing two different types of entangled states. Note 
that the negativity for the two cases can be the same. As a simple example, 
consider the parameter sets I) p = 1/3, 61 = l/-\/2, 62 = 1/v^, ^3 = and 2) 
p = 23/72,61 = 1/v^, 62 = 1/v^, 63 = 1/v^- Each produces a negativity of 
J\f = 5/54. It must also be true for any entanglement measure which provides 
only one number to quantify the entanglement, that there exists parameters 
for which the entanglement is the same, but the types of entanglement are 
different. 

Since the p, and bi are measurable quantities, we may determine the number 
of negative eigenvalues. Alternatively, we could determine number of times 
the coefficients of the characteristic polynomial of p^s change sign. This is 
equal to the number of positive eigenvalues. (See [20,21].) Thus the number 
of negative eigenvalues of the partially transposed density operator can be 
extracted experimentally and provide a sufficient condition for distinguishing 
between types of entangled states. 



5 Conclusions 

DPS are simply described in terms of a pure state component and a polar- 
ization length. Each of these states has a large invariant subspace making it 
tractable to compute in closed form several quantities such as distance metrics 
between states and entanglement between subsystems in a joint depolarized 
state. Such quantities are useful for determining the distinguishability of quan- 
tum states and the nature of quantum correlations that could be used for tasks 
such as entanglement distillation. 

Aside from their simplicity, there is a physical motivation for studying such 
states: namely, a continuous subset of such states corresponds to output states 
from physically allowed depolarization channels. Any completely positive map 
can be driven to a depolarization channel by suitable stochastic unitary oper- 
ations, and the strength of the depolarization is dictated by the magnitude of 
the identity component of the map. In this sense the PDPS correspond to the 
output of a standard form of quantum maps with a pure state input. We have 
described how to experimentally measure the parameters of a DPS by mea- 
suring invariants generated by conditions on the coherence vector describing 
the state. Generically, a measurement of all D such invariants on an arbitrary 
quantum state will allow for a complete reconstruction of the spectrum of the 
state. However, given prior knowledge that the state is a DPS (for example 
by beginning with a pure state, applying an unknown quantum map, and de- 
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polarizing), one can obtain the relevant data by simpler means. Specifically 
by measuring two quantities Tr[p^] and Tr[p^], one obtains the depolarization 
strength. For bipartite systems, measurements of the reduced state spectrum 
then allows for a sufficient measure of entanglement between subsystems via 
the negativity. This requires only Da measurements and is a considerable sim- 
plification versus tomography on the joint state. These measurements can also 
be used to find the number of negative eigenvalues of the partially transposed 
density operator. This number can be used to provide qualitative information 
about the type, as well as amount of entanglement present in the joint state. 
This could, for example, help to distinguish between SU(2) and SU(3) singlet 
states thus providing information about the types of interaction between two 
distant objects. 

We have shown that for bipartite systems with composite dimension D = 
DaDbi the negativity of DPS is identically zero if p < l/(-Dmax(6j6j/) + 1). 
Yet it is also known that the state is separable ifp<l/(-D/2 + l). Do there 
exist bound entangled DPS in between? Verifying the existence of bound en- 
tangled states requires searching in the region of positive partial transpose 
states for states which are not separable. This can be done by constructing 
operators which give witness to separability. Many results have been obtained 
for low rank states [30] , but our case is maximal rank (because of the presence 
of the identity component). Recently, work [31] has shown the existence of op- 
timal separability witnesses for a class of three parameter mixed states. These 
states are bipartite systems with equal dimension composed of the identity 
mixed with three maximally entangled states (locally equivalent to the state 
I $''")). The authors numerically find bound entangled states when two of the 
parameters are nonzero. It is possible that this analysis could also assist in 
finding, or ruling out, bound entangled DPS. 



A Examples 

The examples of this appendix show some reproductions of known results using 
our simplified methods. The first provides a canonical form for DPS for all two 
qubit states. The second example shows that the limits for the separability of 
isotropic states can be derived using our methods given that isotropic states 
subset of the DPS. 

A.l Two Qubits 

This section contains an explicit example of two qubits which are in DPS 
form. The example shows the reduced number of parameters-just two relevant 
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parameters for two qubits determined up to local unitaries-which are present 
in a DPS. 



If a pure state represented by p is acted upon by a depolarizing channel and 
the original po was a singlet state for two qubits, the result of the depolarizing 
channel is called a Werner state [38] . For a four dimensional Hilbert space, the 
density matrix has the form [19,20,21] 

p= i (l + V6n-A) . (A.l) 



This may also be written, for a two-qubit system as 

p = i(l ® 1 + 5^ ■ ® 1 + 1 ® 5^ ■ + a.af ® af ). (A.2) 
4 ■' 



In this expression the constant factor, \/Q has been absorbed into the expan- 
sion coefficients 5*"^, which represent the "spins" of the first and second 
particles respectively and Cjj which represents the correlations between par- 
ticle states. In Ref. [39] an explicit canonical form is given for a two qubit 
pure-state density matrix in the Schmidt form: 



p= ^(1 ® 1 + cos(fi)(T2 ® 1 + cos(r2)l ® Oz 

-|- sin(n)(T^. ® Ox — sin(f2)(Tj^ ® (Jy^ (J^® cr^). 

(A.3) 

where < f2 < 7r/2. The original density matrix is related to this one by a 
set of local unitary transformations. The fact that an explicit canonical form 
has been given for the pure state density matrix of two qubits implies that an 
explicit form of a depolarized pure state density matrix can also be given. 

For a DPS the canonical form is given by 



Pd = -(l®l+p cos(f2)cr2 ® 1 + p cos(n) 1 ® o"^ 

-|-psin(r2)cra; ® — psin(f2)(jy ® (Jy +po'z ® cTz)- 

(A.4) 

The partial transpose (with respect to either subsystem) gives the following 
eigenvalues, 

11 11 ^ > 

= I'^^-P) + 2PSin{n)- p4 = 4(1 -p) - 2psm(fi). 
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The only eigenvalue which could be negative is /i4. However, cannot be 
negative for p < 1/3 = 1/(1^/2 + 1) which is the condition derived earlier. 
Otherwise, ii p> 1/3 we can have entangled states when 

sm{n) > (A.6) 



Note that, if p < 0, the roles of /i2(Ati) and /i4 (/is) are reversed and the same 
conditions apply with the added condition that the original density operator 
is not positive if p < —1/3. 

One might refer to these as generalized Werner states since, unlike the original 
Werner states, there are two variable parameters. One describes the magnitude 
of the coherence vector and the other describes the entanglement of the pure 
state to which the DPS corresponds. 



A. 2 Isotropic States 



We can verify an entanglement property of isotropic states with this result. 
Isotropic states are defined over bipartite states of equal dimension, {Da = 
Db) by a single parameter F 



Pd 



n2 



^AB 



AB 



AB 



where \^^) ~ 75X^^=0 ^ ^-^^^ \3)b- "^^^^ class of states parameterizes de- 
polarized maximally entangled states where p = {D\F — 1)/{D\ — 1). The 
entanglement properties of isotropic states have been studied before and it 
has been shown [40] that pf is separable iff < F < 1/Da, or —1/(D\ — 1) < 
p < 1/{Da + 1). This is consistent with the above result as all pairs of 
Schmidt coefficients for pp have the value bjbj' = I/Da which means that 
for p > 1/ {Da + 1), the isotropic states are entangled. Consequently, there are 
no bound entangled isotropic states. 
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